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The zero-temperature phase diagram of a Bose-Einstein condensate confined in realistic one- 
dimensional ^-periodic optical superlattices is investigated. The system of interacting bosons is 
modeled in terms of a Bose-Hubbard Hamiltonian whose site-dependent local potentials and hopping 
amplitudes reflect the periodicity of the lattice partition in l-site cells. Relying on the exact mapping 
between the hard-core limit of the boson Hamiltonian and the model of spinless noninteracting 
fermions, incompressible insulator domains are shown to exist for rational fillings that are predicted 
to be compressible in the atomic limit. The corresponding boundaries, qualitatively described 
in a multiple-site mean-field approach, are shown to exhibit an unusual loophole shape. A more 
quantitative description of the loophole domain boundaries at half filling for the special case £ — 2 
is supplied in terms of analytic strong-coupling expansions and quantum Monte Carlo simulations. 
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Ultracold bosons trapped in optical lattices provide a 
direct realization of the Bose-Hubbard (BH) model origi- 
nally introduced to describe liquid He in confined geome- 
tries P, • In this framework, the sites of the ambient 
lattice correspond to the local minima of the effective 
potential produced by the interference pattern of coun- 
terpropagating laser beams. The bosons move across the 
lattice with a hopping amplitude whose simple relation 
with the intensity of the laser beams 0] allows an un- 
precedented tuning of the model parameters. Such con- 
trol played a key role in a recent breakthrough experi- 
ment by Greiner and co-workers who were able to 
reveal the superfluid-insulator quantum phase transition 
characterizing the Bose-Hubbard model |jy]. 

The superposition of optical lattices with different peri- 
ods allows to generate periodic trapping potentials char- 
acterized by a richer spatial modulation, the so-called 
optical superlattices 0, 0, S 0> 13 • As it is well known, in 
general the single particle spectrum of superlattice (SL) 
models features band-gaps that critically affect the phys- 
ical properties of the system. In this paper we investi- 
gate the zero-temperature phase diagram of the superlat- 
tice BH model, showing that classes of parameter choices 
exist giving rise to unusual loophole-shaped fractional- 
filling insulator domains. In order to reduce unessential 
details, we focus on the emblematic case of realistic one- 
dimensional superlattices with simple periodicity, that 
can be created following the scheme described in Ref. 4 
In this case, the presence of insulator loophole domains is 
directly related to the band-gaps occurring in the single- 
particle spectrum of the model. The extension of our re- 
sults to generic superlattices will be presented elsewhere. 

The BH Hamiltonian on a ID superlattice comprising 



M sites reads 

= ^ [ -^nk{nk - 1) - (/i - Vk)nk 

k=l 

- tfc(afcafc+i + a^afe+i) (1) 

where a^, ak and Uk — a'^a^ are respectively the boson 
creation, annihilation and number operators relevant to 
the site labeled k. As to the Hamiltonian parameters, 
[/ > accounts for on-site repulsion (proportional to 
the atomic scattering length), /i is the grand canonical 
chemical potential, Vk is the local potential at site k and 
tk is the hopping amplitude between nearest neighbour 
sites k and k -\- 1. According to the tight-binding-like 
approach of Ref. 13, on ^-periodic optical superlattices 

= tT^-nod(kl) + li Vk = l't^mod(fe,£) + l7 (2) 

where Th and (with h = 1, 2, . . . , £) describe the mod- 
ulation of the hopping amplitude and local potential 
within each cell, while t and v are scaling factors. 

The zero-temperature phase diagram of Hamiltonian 
(fT|l can be characterized by the variation in the ground- 
state (free) energy produced by a finite variation of the 
number of bosons in the system. In details, the sys- 
tem is either in an insulator or superfiuid state depend- 
ing on whether such energy variation remains finite or 
vanishes in the thermodynamic limit M — > oo. Note 
that, since the number operator commutes with 

Hamiltonian the ground-state energy of the latter 
can be conveniently studied in the subspaces relevant to 
a fixed number of particles. Denoting E{N) the gromid- 
state energy of H restricted to the A^-boson subspace, a 
critical incompressible domain exists whenever the equa- 
tion E{fM) — E{fM ± 1) yields two distinct solutions 
/i_ ^ in the thermodynamic limit. In this case /i± 
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give the boundary of the insulator domain relevant to the 
critical filling / = N/M in the n/U — t/U phase diagram. 
Note that the compressibility of the system, k — d^j^N, 
vanishes within these domains. 

On simple lattices, Vh = v = 0, — t — 1, the pres- 
ence of the energy gaps characterizing the incompress- 
ible insulator domains can be explained intuitively based 
on the so-called atomic limit, t —>■ 0, where the sites of 
the lattice can be regarded as independent "atoms" 0. 
In such limit it is easy to prove that the free energy is 
always minimized by an integer filling. More precisely 
TV = for ^ < 0, while N = qM for q < n/U < q+l, 
where g G N. Notice that the entire set of non commen- 
surate populations qM < N < {q + 1)M equivalently 
minimizes the free energy for the singular value of the 
chemical potential /i = qll. The gain in kinetic energy 
produced by switching a small hopping amplitude on is 
sufficient to make these configurations favourable on a 
narrow fi interval around qU , and the interval pertaining 
to the integer-filling configuration consequently shrinks 
by a small amount. Several numerical and analytical ap- 
proaches (See Ref. 3 and references therein for a review) 
show that the region of the n/U — t/U enclosed by the 
boundaries /i± has a lobe-like shape whose height de- 
creases with increasing /. 

In general, on ^-periodic superlattices incompressible 
domains occur with critical rational fillings / = N/M — 
q/i, where q & N [ESUl- This is quite easily understood 
in the extreme situation where the cells are drastically 
isolated from each other, i.e. when the hopping ampli- 
tudes between neighbouring_cells are much smaller than 
those within the same cell 7]. This is basically a super- 
atomic limit, the isolated cells being the super-atoms. In 
this situation the rational- filling f = q/£ oi the system 
corresponds to an integer-filling f = q in terms of super- 
atoms. A similar argument applies when the cell features 
uneven local potentials Vk 0, @ ■ Again, in the atomic 
limit t —> 0, the rational filling f = q/£ corresponds to an 
integer-filling f = q for the energetically more favourable 
subset of atoms obtained extracting from each cell the 
mod(g, £) + l sites featuring the lowest local potentials. In 
general, assuming that max({w/j}) — min({ufc}) < U and, 
without loss of generality, that inhi({vk}) = 0, the chem- 
ical potential interval q < (i/U < g -|- 1, g S N, is divided 
into t parts, the j-th subinterval being the basis of the in- 
sulator domain relevant to the rational filling / = q+j/(. 
The points marking the boundaries of these intervals are 
directly related to the local potentials, /i — NU + Vj, 
j = 1, . . . , £ 13]. The condition max({t;fc}) > U, discussed 
to some extent in Ref.j^ makes things more complex, and 
we do not address it here. Note that whenever the local 
potentials Vk are not all different from each other, some 
of the above intervals reduce to singular points, and no 
insulator domain seems to correspond to the relevant ra- 
tional fillings. 

In this paper we show that the zero-temperature phase 
diagram of a BH model on a ID ^-periodic SL features 
critical q/£ rational filling insulator domains for every 




FIG. 1: Zero-temperature phase diagrams for two simple 
£ = 3 superlattices (sketched in the insets). Upper panel: 
Ti — T2 = T3 = 1, vi = 0, V2 = V3 = 0.5, V = U. Lower 
panel: ti = T2 — 1, — 0.3, i^i — U2 = i^s ~ 0. 



q G N. That is to say, an incompressible insulator domain 
exists also for those rational fillings that are predicted to 
be compressible in the atomic limit. We prove this by ex- 
ploiting the exact mapping between the hard-core limit 
{t/U — > 0) of Hamiltonian and free spinless fermions 
in one dimension |lOl |. We furthermore investigate the 
boundary of these domains showing that they feature an 
unusual loophole shape. The presence of loophole insula- 
tor domains is qualitatively predicted by a multiple-site 
mean-field approach generalizing the one introduced in 
Ref. [m Quantitative results are obtained via quantum 
Monte Carlo (QMC) simulations. Furthermore, an ana- 
lytic perturbative description of the boundaries for the 
half- filling insulator domain of the special case £ = 2, 
vi = V2 = is obtained which proves quite satisfactory. 

As we mention, the hard-core limit, t/U <^ 1, 
of ID Hamiltonian ^ can be mapped exactly onto 
the tight-binding Hamiltonian describing spinless non- 
interacting fermions on the same ambient lattice |lClj| . 

= - Vk)c'lck - tk(ckc'/:^^ -I- Cfc+ic^)], where 

Ck is the fermionic annihilation operator and {cfc,c^} = 
5kh E2- Since the fermions do not interact, the zero- 
temperature many-body ground state of iJp is simply 
obtained in terms of corresponding single particle spec- 
trum. As it is well-known, on a ^-periodic SL, the lat- 
ter features £ disjoint bands, each containing M /£ states 
[l3j |. Hence, for critical rational fillings q/£ {q — 1, . . . ,£) 
the fermionic system features an insulating ground-state 
which corresponds to an incompressible Mott state of the 
bosonic system. This proves that an insulating domain 
exists for any rational filling, including the ones that are 
not critical in the atomic limit. In particular, the latter 
feature a cusp rather than a wide basis ioi t/U 0. 

Qualitative information about the shape of these do- 
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FIG. 2: Left and right panel show a magnification of the first 
loophole insulator domains appearing in the upper and lower 
panels of Fig. respectively. 



mains can be gained by considering a multiple-site mean- 
field approximation jlil to Hamiltonian . We recall 
that, following Ref. Illl the latter can be recast as the 
sum of M coupled single-site Hamiltonians by assuming 
that, for every fc, 



(afc>afc+i + «fc(afe+i) - 2(afe)(afc+i), (3) 



where (•) denotes expectation value on the ground state. 
This approach gives a qualitative description of the frac- 
tional insulator lobes when the local potentials within 
a cell are different from each other jsj, but fails to pre- 
dict the existence of loophole insulator domains. If the 
non-trivial periodicity of the SL is taken into account 
by adopting approximation (j2Jl every £-th site, Hamilto- 
nian can be recast as the sum of ^-site Hamiltoni- 
ans. Note that the ensuing self-consistency problem is 
much more demanding than the usual single-site mean- 



field approach. Indeed, in general it requires the iterative 
diagonalization of a matrix whose size is X]^=o ^q](i Zi ji 
where Q provides a cutoff for the infinite Hilbert space of 
the problem |[l5j , and must be sufficiently larger than the 
typical population for a given choice of the Hamiltonian 
parameters [T^ . We mention that a two-site mean- field 
approach is adopted in Ref. for the study of homoge- 
neous lattices. 

Figure ^ shows zero-temperature phase-diagrams for 
two simple £ = 3 superlattices featuring insulator loop- 
hole domains, as evaluated by the multiple-site mean- 
field self-consistent approach. As usual, the insulator do- 
mains (white areas) are defined by the vanishing of the 
superfluid order parameter, (uk) at every site k of the SL. 
The relevant effective optical potentials are schematically 
depicted in the insets. Figure [3 shows a magnification of 
the loophole insulator domains appearing in Fig. ^ 

A more precise approximation of the insulator do- 
main borders can be achieved by means of perturba- 
tive approaches. On homogeneous lattices E{fM) and 
E{fM ± I) can be evaluated perturbatively, t being the 
perturbative parameter Il7l|. A reasonable effort allows 
to extend this approach to superlattices, for the insula- 
tor domains that are predicted by the atomic limit i — > 0. 
Analytic results for the entire class of superlattices with 
constant local potential, where the atomic limit predicts 
integer filling domains only, are reported in Ref. [13 ■ 

Similarly to the multiple-site mean-field approach, the 
loophole domains can be perturbatively studied by taking 
into account the periodicity of the SL, i.e. assuming that 
the perturbative parameter is the hopping term between 
adjacent cells. More precisely, if min({Tfc}) — Tg, one can 
write 



M Ml I M/l-H-i 

k=l s=l s=0 h=l 



and consider ti as the perturbative parameter. The re- 
sults are expected to be satisfactory for SLs very close to 
the super-atomic limit, i.e. when T£ ^ Tfc, V/c ^ t. This 
approach can be carried out analytically for the lowest 
loophole domains of some simple SLs. For instance, the 
first-order analytical result for the half-filling loophole 
insulator domain in the simple case £ — 2, vi — V2 'is 



{T2-Tl)t 



T2 {AtTi +U + RY 

2 16(iTi)2 ^u'^ + UR 



U-R 



(5) 
(6) 



where R = sJU"^ + X^itT^Y . Fig. O shows the loophole 
domain borders as obtained by Eqs. and © for the 



parameter choices ri = 1 and = 0.2 (solid curves) and 
T2 = 0.6 (dashed curves). The data points in Fig.|31are 
the numerical QMC estimates of fi± for a SL comprising 
M — 100 sites. First of all we notice that these data 
points quantitatively confirm the loophole shape of the 
fractional incompressible domain under concern. On the 
other hand they show the perturbative results to be quite 
satisfactory also far from the super- atomic limit, T2 <C ti. 

The algorithm we adopted for the numerical evaluation 
of the data points in Fig. O is the so-called population 
QMC, which basically amounts to a stochastic version of 
the power method for finding the maximal eigenvalue of 
a matrix p^ . 

In summary, in this paper we prove that the zero- 
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FIG. 3: Inset A.: Schematic representation of the optical 
potential characterizing the ^ = 2 SL under concern. The rel- 
evant parameters are i/i = 1^2 = {v = 0), ti ^ T2. Inset B.: 
Boundary of the half-filling loophole insulator for the param- 
eter choice Ti = 1, tau2 — 0.6. Solid curves: analytical result 
according to first order perturbative expansions, Esq. Q and 
®; Data points: QMC simulations. Main plot: Boundary 
of the half-filling loophole insulator for the parameter choice 
Ti = 1, T2 = 0.2 according to the perturbative expansion 
(dashed curves) and QMC simulations (data points). The 
perturbative result of inset B. is also shown for comparison 
(solid curves). The broadening of the QMC errorbars with 
increasing t/U corresponds to the increase in the number of 
Fock states contributing to the ground state, which makes the 
estimate of the energy more noisy 19]. Furthermore, finite- 
size efi'ects introduce an upper bound on the precision of the 
estimate at the tip of the domain. 



temperature phase-diagram of a ^-periodic ID superlat- 
tice BH Hamiltonian features a g/f- filling incompressible 
insulator domain for every q € N. Furthermore we show 
that the insulator domains that are not predicted in the 
atomic limit feature an unusual loophole shape. It is 
worth remarking that such shape implies quite interest- 
ing reentrant behaviours. As it is clear from Fig. |21 for 
some values of the chemical potential, the system rather 
unusually enters an insulating phase as a result of an 



increase in kinetic energy. A similar behaviour is ob- 
served also in the first lobe of the ID simple lattice (2(1] . 
Even more unusually, two distinct insulator domain sepa- 
rated by a superfluid stretch can occur at the same /i/C/. 
The presence of the loophole insulator domains can be 
possibly detected as described in Ref. 3]. Recently, an 
alternative scheme has been proposed, based on the mea- 
sure of the Mott domain width via the laser probing of 
the particle-hole energy gap characterizing the insulat- 
ing phase [23|- This scheme seems suitable for detecting 
the reentrant character of the loophole domains. Indeed, 
as a function of t/U, the measured width is expected to 
feature a maximum for loophole domains, whereas it is 
simply non increasing for the usual lobes [2ll ]. 

We emphasize that the above considered ID superlat- 
tices are quite realistic. Indeed two superposing laser 
beams give rise to a ID optical lattice whose period is 
determined by their crossing angle When two such 
lattices with commensurate period are superimposed, a 
ID optical superlattice is obtained Q. This effective 
trapping potential can be used to fragment an elongated, 
cigar-shaped Bose-Einstein condensate, giving rise to a 
system described by Hamiltonian Q • Ref. Q reports on 
the realization of a € = 3 ID optical superlattice where 
ultracold bosons are trapped. Note furthermore that the 
optical potentials in Fig. ^ basically differ by an overall 
sign. This means that in principle it is possible to switch 
between them by changing the sign of laser detuning. 
We mention that the phase diagrams in Figs. CJ-© are 
obtained assuming that only t is varied, for instance by 
varying the laser beam intensity /, while the remaining 
parameters in Eq. (0) remain fixed. This is quite natural 
when there is no spatial modulation in the local poten- 
tials, and one can safely set = 0. This not being the 
case, the linear dependence of w on / can be ignored to 
a first approximation, since t ~ e . If this is not sat- 
isfactory, the desired effect can be obtained by adjusting 
both the laser intensity and the beam setup. 
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